Higher correlations, universal distributions and finite size scaling in the 
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Recently we constructed a renormalizable field theory up to two loops for the quasi-static depinning of elastic 
manifolds in a disordered environment. Here we explore further properties of the theory. We show how higher 
correlation functions of the displacement field can be computed. Drastic simplifications occur, unveiling much 
simpler diagrammatic rules than anticipated. This is applied to the universal scaled width-distribution. The 
expansion in d = 4 — e predicts that the scaled distribution coincides to the lowest orders with the one for 
a Gaussian theory with propagator G{q) = l/g'*^^^, ^ being the roughness exponent. The deviations from 
this Gaussian result are small and involve higher correlation functions, which are computed here for different 
boundary conditions. Other universal quantities are defined and evaluated: We perform a general analysis of the 
stability of the fixed point. We find that the correction-to-scaling exponent is uj = —e and not — e/3 as used 
in the analysis of some simulations. A more detailed study of the upper critical dimension is given, where the 
roughness of interfaces grows as a power of a logarithm instead of a pure power. 



I. INTRODUCTION 

Understanding the behavior of an elastic interface in a ran- 
dom potential is important for many experimental systems and 
still offers a considerable theoretical challenge [1^]. It is ex- 
pected that below the upper critical dimension due the inter- 
face is pinned by arbitrarily weak disorder, into some rough 
configurations and that at zero temperature it can acquire a 
non-zero velocity under an applied force / only if / is larger 
than the depinning threshold fc- A functional renormalization 
group (FRG) method predicts that due — 4 for the statics [5], 
and for the simplest universality class, the so called isotropic 
depinning [6-8]. 

There has been recent progress towards a precise descrip- 
tion of the depinning transition. From the theory side, the 
FRG for single component manifolds, originally studied to 
one loop in an expansion in e = due — d, has now been ex- 
tended to a field theory shown to be renormalizable to two 
loops. Renormalizable, we recall, means it has a well defined 
continuum limit, which is independent of all microscopic de- 
tails, and thus ensures universality of large scale observables. 
Presumably there exists a fully renormalizable theory to all 
orders, with full predictive power [9-11]. From the side of 
numerics a novel powerful algorithm allows to obtain the con- 
figurations at (or just below) depinning with much improved 
accuracy [12-14]. A reasonable agreement between the two 
methods was found in a measurement of the roughness expo- 
nent (, especially the clear conclusion that ( > e/3 contrarily 
to a previous conjecture [7, 8] (C = e/3) based on the 1-loop 
analysis. 

The field theory of depinning in its present form is uncon- 
ventional in that one must work with a non-analytic action. 
This peculiar feature is a deep part of the physics of the prob- 
lem and necessary to avoid the so called dimensional reduc- 
tion. It makes the perturbation theory superficially "ambigu- 
ous". A non-trivial step taken in [9-11] to define the theory 
at T = as the limit u — > of the moving phase, was to 



assume that the interface-position is monotonic in time. This 
removes the ambiguity and, remarkably, leads to a renormaliz- 
able theory, to at least two loops [9-11]. This is supported by 
the "non-crossing theorems" which apply to single component 
depinning and, remarkably, is the same property allowing to 
show ergodicity and to construct an efficient algorithm to find 
the exact critical configuration at depinning [12, 15]. The ori- 
gin of recent progresses in both numerics and field theory are 
thus related. Clearly one would like to test this novel field the- 
ory by calculating more universal measurable quantities and 
study its properties. 

In this paper we further explore the field theory constructed 
in [9-11]. We study displacement correlations of more than 
two points. We find that these correlations are static. Al- 
though physically natural, if one wants quasi-static depinning 
to make sense, this manifests itself through rather non-trivial 
massive cancellations in the time dependence of multi-point 
diagrams. We elucidate these cancellations and obtain as a 
consequence for a large class of diagrams much simpler di- 
agrammatic rules than previously anticipated. Basically, all 
time integrals become almost trivial, resulting in a theory with 
"quasi static" diagrams. We then apply these properties to the 
calculation of universal observables. One natural universal 
quantity is the so-called width distribution of the interface. 
Interestingly, to the two lowest leading orders in e — 4 — d, 
the distribution coincides with the one for a Gaussian theory 
with the full non-trivial propagator G{q) = 1/q'^^'^'^, C be- 
ing the depinning exponent. This is also the subject of a re- 
lated publication [16], where the distribution is also measured 
numerically. Here we give a detailed presentation and also 
compute the higher connected cumulants of the displacement 
field, i.e. deviations from the Gaussian. Some of these results 
are quoted in [16]. 

In a second part we study the theory at the upper critical 
dimension. The motivation is that no exact result is avail- 
able to confirm that due = 4 (the only exactly solved limit 
corresponding to fully connected models [17-19]). Thus the 
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question of what is the upper critical dimension d^c is still de- 
bated, even though the field theory of depinning [6-11, 20] 
clearly predicts due = 4. Also, in the other class of depin- 
ning transitions, the so called anisotropic depinning class with 
KPZ nonlinearities, there is not even a convincing prediction 
for due [21-23] and recent numerical studies have reopened 
the debate [24]. Recently it has become possible to study nu- 
merically depinning and statics in high dimensional spaces for 
reasonable system sizes with better precision, allowing for the 
hope to settle the issue of the upper critical dimension in the 
near future [12-14, 25]. It is thus important to give precise 
predictions for the behavior predicted by the FRG, in order to 
compare with numerics. 

Finally, we also clarify the issue of finite size scaling. In 
a previous work, used in several simulations, the value w = 
— e/3 was used for the finite size scaling exponent [26, 27]. 
We find that the correct value is instead w = —e. This may 
prove useful in numerical studies [35]. 

The paper is organized as follows. In Section 11 we define 
the model, briefly review the FRG method and field theory 
and define the main observable of interest here, the width dis- 
tribution. In Section III we compute the Laplace transform 
of the width distribution in perturbation theory and find that 
to lowest order in e it coincides with a Gaussian Approxima- 
tion. This approximation is introduced and further studied. 
Some results on Laplace Inversion are given. In Section IV we 
go beyond the Gaussian Approximation and compute higher 
connected cumulants. The detailed calculation of the fourth 
cumulant (4-point connected correlation function of the dis- 
placement field) is given, and the cancellations that occur in 
the field theory are studied. In Section VI we discuss the up- 
per critical dimension and in Section VII the finite size scal- 
ing. The effect of various boundary conditions is studied in 
Appendix A. 



II. MODEL AND OBSERVABLES 
A. Model 

We study the over-damped dynamics described by the equa- 
tion of motion 



vdtUxt = cVluxt + F{x, Ua;t) + f 



(2.1) 



with friction rj. Long range elasticity relevant for solid friction 
at the upper critical dimension d = 2, can be studied replac- 
ing cq"^ c\q\. In presence of an applied force / the center 
of mass velocity is f = L~'^ dtUxt- The pinning force is 
x) = —duV{u, x) and thus the second cumulant of the 
force is 

F(x, u)F{x', u') = A{u - u')6'^{x - x') , (2.2) 

such that A(u) = —R"{u) in the bare model, where R{u) is 
the correlator of the random potential. Random bond disor- 
der is modeled by a short range function R{u), random field 
(RF) disorder of amplitude a by R{v) ^ — cjul at large u and 
CDW disorder by a periodic function R{u). 



B. Review of FRG and field theory 

Let us briefly review the field theoretic approach, more details 
can be found in [10]. The dynamical action (MSR) averaged 
over disorder is given by e"*^ with 



S[u,u]= Uxt{dt-yDuxt-7i I UxtA{Uxt-Uxt')Uxt' 
Jxt ^ Jxtt' 

(2.3) 

Here and below we denote ■= J d'^x, in Fourier J^. := 
/ and = / dt. The FRG shows that the full func- 
tion A(u) becomes relevant below d = due = 4 and a flow 
equation for its scale dependence has been derived to one and 
two loops, in an expansion in d = 4 — e. In Ref. [10] this 
was derived by adding a small mass m as an infrared cutoff 
and computing the flow of disorder, defined from the effective 
action r[u, ii] of the theory, as m decreases towards zero. As 
in [10] we will denote by Ao(m) the bare disorder correlator, 
i.e. the one appearing in the action S in (2.3), and by A(u) 
the renormaUzed one, appearing in T[u, u] which has a simi- 
lar expression as (2.3). The rescaled disorder is then defined 
by 



A(u) = ^rn^"^^A(um'^) , 
eh 



(2.4) 



where Ii 



= + m^)^^ is the 1-loop integral. 

It was then shown in [9, 10] that (2.3) leads to a functional 
renormalization group equation 

-m9„A(u) = (e - 2C)A(«) + CuK'{u) 
1 



{Mu)-m)\ 

{K{u)-m)K\uf 



(2.5) 



up to O(A^) terms. This equation implies that there are only 
two main universality classes at depinning, a single RF fixed 
point for interfaces and a periodic one for CDW type disor- 
der [9, 10]. Both ^ and the fixed point function A*{u) were 
determined to order O(e^) for these classes [9, 10]. 

The important feature of the field theory of depinning is that 
A(?i) has a cusp-like non-analyticity at u = 0. As was shown 
in [9, 10] calculations in the non-analytic theory (e.g. yielding 
(2.5)) are meaningfully performed using the expansion: 



A(?i) = A(0) + A(0+)|u| + ^A{0+)u^ 



(2.6) 



Performing Wick averages yields the usual diagrams, except 
that their actual values involve averages of e.g. sign func- 
tions of the fields. Replacing everywhere sgn(ut — Uf) 
sgn{t — t') is justified for single component quasi-static de- 
pinning (i.e. in the limit of vanishing velocity v = 0+). This 
yields diagrams with sometimes complicated internal time 
and momentum dependences. We find however that in some 
cases massive cancellations occur despite the complications 
due to the time dependence between various diagrams, con- 
tributing to the same observable. 



3 



C. Universal distributions and observables 

To motivate the present study let us consider one specific ex- 
ample of a universal observable, the width distribution of the 
configuration at depiiming (the so-called critical configura- 
tion). The width of a configuration is defined in a given disor- 
der realization as 



w 



uf 



(2.7) 



where u= J^Ux is the center of mass and L"^ the volume 
of the system. The basic observation is that the sample to 
sample probability distribution P{w'^) of is expected to 
be universal, with a single scale set by the disorder averaged 
second cumulant w"^, i.e: 



1 



P{w') = =f = 



w 



(2.8) 



f{z) is a universal function. This holds for thermal averages 
in a number of finite temperature problems of pure systems 
[28, 29]. Here we show that it also holds for depinning at 
T = and compute the distribution, first within a simple 
Gaussian approximation and second within the e-expansion. 
In the process we study higher point correlation functions in 
the field theory of depinning, define specific universal ratios 
of these, describing deviations from Gaussian behavior and 
compute them. 

Before turning to actual calculations let us first summarize 
the general spirit of the method and discuss the question of 
the universality of such observables. The hallmark of a renor- 
malizable theory is that if one expresses the correlation func- 
tions in an expansion in the renormalized disorder A, then 
the resulting expressions are UV finite, equivalently they have 
a well defined continuum limit, independent of short scale 
details. On a technical level, this can be achieved by com- 
puting correlations in standard perturbation theory to a given 
order in powers of Aq, and then using the relation between 
renormalized disorder A and bare one Aq to the same order, 
or equivalently through the definition of appropriate counter- 
terms [36]. Here, we restrict ourselves to calculations at dom- 
inant order in e and thus using either A or Aq makes no dif- 
ference. Beyond the Larkin scale, however, these are non- 
analytic functions, which is crucial. 

In the limit of large scales or large system sizes, the fixed 
point form reached by the rescaled A implies that the resulting 
observable, e.g. the width distribution, is universal. Universal 
means that these quantities do not depend on the short scale 
details. However they do depend on the details of the large 
scale infrared (IR) cutoff, i.e. of the type of chosen boundary 
conditions. Here we focus on periodic boundary conditions 
[34], also of interest for numerical simulations [16], although 
we sometimes give results for other types, for instance for the 
massive IR cutoff described in the previous paragraph. 

Since the FRG method developed in [9, 10] and summa- 
rized above uses a mass as IR-cutoff and defines disorder ver- 
tices at zero momentum, one should be careful in calculations 
with e.g. periodic boundary conditions. Since we only com- 



pute observables either to dominant order in e, or within a 1- 
loop approximation, it is easy to make the necessary replace- 
ments, as will be indicated below. For instance the 1-loop 
FRG equation remains identical to the two first fines of (2.5), 
the only changes being that 

1 . — m9m A has to be replaced by A. 

2. m ^ 1/ L'm the definition of the rescaled disorder 



3. the 1-loop integral I\ = 



^ entering into the 



k (fe^+m^) 

definition of the rescaled disorder has to be replaced by 
its homologue for periodic boundary conditions [37]: 



used below. 



(2.9) 



III. WIDTH DISTRIBUTION: 
PERTURBATION THEORY AND GAUSSIAN 
APPROXIMATION 

Let us start by giving the simplest approximation for this dis- 
tribution. It can be derived in two ways: (i) perturbation the- 
ory in the renormalized theory to lowest order in e (ii) a simple 
Gaussian approximation. In the end this will motivate going 
further, i.e. studying deviations from the Gaussian approxima- 
tion. 

A. Perturbation tlieory 

We now study perturbation theory. To compute the width dis- 
tribution using the dynamical field theoretic method [9, 10] 
one can start from the Laplace transform 



W{\) 



(3.1) 



with w"^ = Ylixi^x — uf. Here and below we omit the global 
multiplicative factor L~'^ in the definition of w"^, since in the 
end we will always normalize the distribution P{w^) by fix- 
ing its first moment to unity (in (3.1) it can be absorbed by a 
rescaling of A). 

Expanding in powers of the correlator of the pirming force 
A(u) (to lowest order this is equivalent to Aq see above), one 
finds that to leading order In W{X) is the sum of all connected 
1-loop diagrams, as represented in Fig. 1. The loop with n 
disorder vertices and n insertions of w"^ is 



2n 



E 



-2AA(0) \" 

{q^r ) 



(3.2) 



Here and below the sums over q thus runs over a d- 
dimensional hypercubic lattice with spacing and the 0- 
mode is excluded, as appropriate for periodic boundary con- 
ditions. If one were now to resum (3.2) over n one would 
find: 



W{X) = l[il + 2XGiq)y 



-1/2 



(3.3) 
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r.h.s. of (3.5) generalizes into 



— □ — e^f 





FIG. 1: Examples of contributions to (3.1), terms , (w^)"^ and 
(w^)^ (bottom), together with the vertices for disorder (top left), 
(top center), and response function (top right). 



with G{q) = GLarkinC?) = A(0)/g^, a Larkin model type 
result if interpreted as naive perturbation theory (i.e. if A(0) 
was interpreted as the bare original disorder rather than the 
renormalized one). The correct procedure implies that A(0) 
is the running renormalized disorder A(0) A((0) — 
(e/i)-ie(2«-<^)'A*(0) where A*(0) is the (non-universal) 
value of the fixed point [9, 10]. Although I = ln(L) for the 
zero momentum disorder vertex, one notes that a momentum q 
flows in each vertex and one should take care of this by setting 
I ln(l/g). This yields finally (3.3) with 



Gig) = C/q- 



(3.4) 



where the value of C is non-universal and fixed by w"^ . As wiU 
become clear in the following section, the appropriate choice 
for G{q) is the 2-point finite-size scaling-function G^iq) = 
C/q''-+'^ig[qL) with g{0) = 1. The difference between the 
two above-mentioned choices for / simply amounts to the two 
different limits of small, or large qL. However to lowest order 
in e = 4— d they are identical. (For a calculation of the scahng 
function to next order in e see Appendix J in Ref. [10].) 

B. Gaussian Approximation and beyond 

A more general approach consistent with the previous calcu- 
lation is the following. We first note that the above result (3.3) 
would be exact if the distribution of the displacement fields u 
were Gaussian. It can thus be caUed the Gaussian Approxi- 
mation (GA). To understand why it was obtained here let us 
consider simply the second connected cumulant of the WD 

c 2 

(w^)^ = — (u>2) . This cumulant however is not con- 

nected w.r.t. the u, and thus there is an exact relation: 



(w;2)2 



Jxy 



(3.5) 



where here Gxy = u^Uy is the exact disorder averaged 2-point 
function. The first term is just Wick's theorem and would be 
the full result if the measure of the u were Gaussian. Analo- 
gous formulae exist for higher cumulants: the first term on the 



IGA 



2"-i(n- 1 



./xi,..a 



GxiX2^X2X3 



. . Gx 



(3.6) 

as a simple consequence of Wick's theorem. This is again 
easily resummed into (3.3) which would thus be exact if the 
measure of u were exactly Gaussian. Note that all results of 
[28, 29] for pure Gaussian theories can also easily be obtained 
by the present resummation method (temperature replacing 
disorder). Thus in the GA, the G{q) appearing in (3.3) is the 
exact 2-point function. It can be tested in a simulation by in- 
serting the measured 2-point function in (3.3). In the large 
but finite system size hmit it takes the scahng function form 
Gl((j) discussed above. 

When compiulng to the numerical results for the width dis- 
tribution, it turns out that the GA is a surprisingly good ap- 
proximation even down to c? = 1. This is discussed in details 
in [16]. However, we do not expect the GA to be exact. It is 
thus interesting to compute the deviation D = j'^^ ^^ w^m^ for 

the second cumulant of in (3.5). It is computed below and 
found to be of order while the GA contribution (first term 
in (3.5)) is of order (since G ^ e). Similarly the GA contri- 
bution to (3.5) is 0(e"), while the deviations are found to be 
0(e^"). This can be summarized asu — ^/euo+eui where uq 
is a Gaussian random variable of 0(1) and ui a non-Gaussian 
oneof 0(1). 

Computing deviations from the GA is thus one motivation 
to compute higher point correlations. 

C. Laplace inversion 

Before doing so, let us discuss how the distribution P{w^) is 
obtained through an inverse Laplace-transform, as 



(3.7) 



Noting that in d = 1 (3.3) can also be written as 

W{X) = l[{l + 2XG{q))-\ (3.8) 

9>0 

this is equivalent to 



P{u 



p>0 



Q>0,q¥=P 



G{p)J 



(3.9) 

This formula shows that for large w^, the distribution is dom- 
inated by the first term p = 1, and in practice summing the 
first few terms gives an excellent approximation. 

It is instructive to apply (3.9) to a random walk, where 
G{q) = Xjq^. Using that (n G N) 

^ sin(7rx) 



n>0 



n 1-^. 



one finds in terms of the width 



(3.10) 



(3.11) 



n>0 
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For d> 1 the situation is more complicated. Writing 

P{w^) = / ^ e-'^ n (1 + 2AG(g))-i . (3.12) 

we have e.g. at least multiphcity 2d for each factor in (3.3), 
as long as no component vanishes, but this multiphcity may 
even be higher, as can be seen from the following solutions of 
the diophantic equation (for 2 dimensions) 1^ + 7^ =5^ + 5^, 
6^ + 7^ = 9^ + 2^, and many more. Let us define the class 
C{q) as 

peCiq) if p^^q^. (3.13) 
Let us index these classes by a, and introduce an order 

Ca<Ca' if 9 e Ca and q' € Ca' \q\ < \q'\ (3.14) 



The number of elements of each class is defined as 

\Ca I := number of elements in Ca (3.15) 

We further define 

qa ■= any element out of Ca ■ (3.16) 

Note that since forp e Ca and p ^ (by definition we exclude 
p = 0) \Ca\ is always even. (3.12) can then be rewritten as 

P{w') = / ^ e-^^ 11(1 + 2AG(g<,))-|C«l/2 . (3.17) 

a 

There are poles at A = —[2G{qa)]~^ ■ The sum over these 
poles can after partial integration be written as 



(|C„|/2-l)! \2G{qa) 



|Ca|/2 



d_ 



|C„|/2-l 



H (l + 2AG(g„))- 



\Cc\/2 



(3.18) 



A=-[2G(gc)]- 



IV. HIGHER POINT CORRELATIONS IN 
DEPINNING FIELD THEORY 

In this section we analyze how one can compute higher cor- 
relations in the depinning field theory and obtain simple dia- 
grammatic rules for doing so. These are illustrated on the 4- 
point function. Specific calculations and results will be given 
in the next Section. 

A. Preliminaries 

We want to compute at T = 0, using the dynamical action S 
in (2.3) , e.g. the 4-point expectation value, connected w.r.t. 
disorder (and u) as defined in the previous sections: 

/ ultult' = / (ultult)^ . (4.1) 

Jxy ^ xy 

Similar formulas hold for higher correlation functions. This is 
identical to a connected expectation value w.r.t. the action S, 
denoted hereafter {■ ■ 

The first step is to show that correlations can all be ex- 
pressed as: 

i.'^xt'^yt) ~ I '^xt,x\t\T^xt,X2t2'^yt,xst2,'^yt,Xiti 

Jxi,ti<t 

2:2*2, 2:3^3, 2:4*4) (4.2) 

Here TZ is the exact response function and T the exact effec- 
tive action (sum of IPI graphs) (with the choice e~'" for the 
probability and F^'') is symmetric). This is the standard rela- 
tion between connected correlation functions and the effective 



action (i.e. 1 -particle irreducible vertex functions IVF). The 
simplification here is that a priori the exact 2-point correlation 

function and vertices such as ri^L,-, could also contribute, but 
their contribution vanishes for T = at the depinning thresh- 
old. This is because {uxiUxm) time independent there, and 
then statistical tilt symmetry implies that all IVF's with at least 
one external u leg carrying frequency u> vanish when this fre- 
quency is set to zero (see section II A in [10]). The above 
formula (4.2) generaUzes straightforwardly to any cormected 
2n-point correlation function of the field u in terms of rl^"]j. 

Next one can compute F^'^uu ir* perturbation, using the di- 
agrammatic rules for the non-analytic action arising from the 
expansion (2.6). Let us denote Ea the number of u external 
legs, riy the number of (unsplitted) vertices Aq, nj the num- 
ber of internal lines (response functions) and L the number 
of (momentum) loops. Then one has 2ni, — n/ — Eu and 
L = \ + ni — riy. Here one has Ey, = 4 and thus the lowest 
order contribution has n„ = 4, nj = 4, i.e. it is the 1-loop 
square (since two vertices n„ = 2 imphes n/ = and is 
disconnected, three vertices n„ = 3 imphes n/ = 2 and is 
1 -particle reducible). The two loop corrections are diagrams 
with 5 vertices and so on. Similarly F?"! to lowest order is 
the 1-loop 2n polygon diagram. 

Since two u fields must come out (arrows) of each vertex, 
there are four possible diagrams corresponding to the 1-loop 

square, shown on figure 2. Each line entering a vertex cor- 
responds to one derivative of the A{u) function of the vertex. 
Thus from (2.6) we see that (a) is proportional to A'(0+)'', (b) 
and (c) to A"(0+)A'(0+)2 A(0) and (d) to A"{0+ fA{(]f. 
However, as we will show below using the so-called mount- 



ing property only (a) is non-zero. 

We found two ways to compute diagram (a) (as well as any 
other similar diagrams), a systematic but complicated way, 
and a simple way which uses a very important property and 
not yet fully elucidated of the field theory of depinning, the 
"quasistatic property" described below. To appreciate the ex- 
tent of the cancellations involved in this drastic simplification, 
we start by sketching the systematic method. 

To perform an actual calculation, since each A vertex in- 
volves fields at two times, one must switch to the splitted dia- 
grammatics, as described in [10]. Diagram (a) then becomes 
the sum of 16 splitted diagrams (two choices per vertex) repre- 
sented in Fig. 3 . Note that the last one is zero since it involves 
an acausal loop. That leaves 15 non-zero and non-trivial dia- 
grams. 

These diagrams correspond to the following. One first ex- 
pands <S^/4! using (2.6), which gives, schematically: 

A'(0+)'*. . 

— ^j^j — W1W2S12M12M3U4S34U34M5W6S56W56W7W8S78M78 • 

(4.3) 

In shorthand notations ui = iixi.ti, ^2 = u^-i.i-j' '^12 = ui — 
U2, S12 = sgn(fi — ^2) omitting all space and time integrals. 
On then carries the Wick contractions, yielding 



4! 



2U1U3U5U7S12S34S56S78 



X (-R32 — R42){R54 — Re4){R76 — R86){Rl8 — R28) ■ 

(4.4) 

Rij = (uiUj) is the free response function. The factor of 
2 comes from the two possible time orientations of the loop. 
Expanding the product of response functions yields the 16 di- 
agrams represented in Fig. 3, where space and time labels are 
ordered turning clockwise around the momentum loop. For il- 
lustration let us indicate the full expression of the first diagram 
in Fig. 3, in momentum space: 

^ia.l)iPl2,ti;P23,t3;P3i,h;P41,t7) = (4.5) 

A'(0+)* / sgn(ti - t2)sgn(<3 - t4)sgn{t5 - h) 

J t2t4^tQtg 

XSgn(t7 — ts)Rpj^^t^-t2Rp2,t5-U^P3,t7-t6^PA,ti-t8 



(a) 



(c) 



(b) 



(d) 



FIG. 2: The four 1-loop diagrams with unsplitted vertices which con- 
tribute to the 4-point irreducible vertex function F^*]^^ 




FIG. 3: The 16 1-loop diagrams with splitted vertices which cor- 
respond to diagram (a) in Fig. 2. The last one, which contains an 
acausal loop and thus vanishes, is added here for future convenience. 



with pij = p^ — pj the entering momenta and Rp^r = 
9{T)e^P the free response function. Because of the sign 
functions the evaluation of these integrals, and of the other 14 
non-vanishing diagrams is very tedious and was handled us- 
ing Mathematica. Adding all diagrams, massive cancellations 
occur The final result is very simple and given below. 

Let us now explain the simple method and the properties of 
the theory which lead to it. 

B. Theorems and other properties 

The simple way to compute the 4-point correlations (and 
higher ones) at depinning is based on the very important fol- 
lowing conjectured property: 



(4.6) Quasistatic property 1 : 



All correlation functions 



computed us- 



ing the diagrammatic rules of the quasistatic field theory of 
depinning at zero temperature and exactly at threshold, are 
independent of all time arguments. 

Using relations such as (4.2) for arbitrary times shows that 
an equivalent way to state this property is the following: 

Quasistatic property 2 : 

All T^?"'\{xiti, . . . ,X2nt2n) are independent of 

tl, ■ ■ ■ , t2n- 

This property, which appears as a physical requirement for 
the correct field theory of depinning, implies non-trivial prop- 
erties of the diagrammatics. Although we will not attempt to 
prove it here in full generality, we have checked it on many 
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examples, and believe that it works. We encourage the reader 
to contribute a valid proof. We will however state and prove 
some easier and useful properties below. 

Once the properties 1 and 2 are accepted, evaluation of the 
diagrams drastically simplifies, thanks to the following trick: 
Since the result does not depend on external times, one can 
take these times mutually infinitely separated, with some fixed 
(and arbitrary) ordering. Then one can integrate easily over all 
internal times since the order at each vertex is then specified 
and each sign function has a fixed value. One recalls that in 
the splitted diagrannmatics all non-vanishing T = diagrams 
are trees (see section 11 A in [10]). This can be seen on the 
fifteen non-vanishing diagrams of Fig. 3. Thus integrating in- 
dependently along each tree starting from the leaves yields 
one correlation function per link, since J^, Rq,t-t' = l/?^- 
Performing this calculation on all fifteen diagrams of Fig. 3 
shows that they cancel pairwise since they differ only by a 
global sign, with the exception of the acausal diagram which 
is zero. Thus the final result is the same as if one had kept 
only (—1) times the acausal graph! 

Before giving the final result below, let us now state the 
easier properties. 

Theorem 1 (mounting trick): 

A diagram which contributes to an u-independent vertex 
function is if it contains a vertex, into which no response- 
function enters. 

Examples are diagrams (b), (c) and (d) on figure 2. This 
theorem thus ensures that only (a) is non- vanishing. 

Proof; 

The following figure demonstrates the principle. Note that 
it may be part of a larger diagram. Especially, there may be 
more response- functions entering into the upper disorder. The 
statement is that 




(4.7) 



Since no response-function enters into the lower disorder 
A{u — u'), due to the assumptions this gives A(0), with no 
dependence on time. Thus one can freely integrate over the 
response-function starting at time t. This integral for both di- 
agrams is / dtR{k,t) = The difference in sign comes 
from deriving the two different ends of the upper disorder ver- 
tex. Thus both contributions exactly cancel. 

Thus T^{^lf^,^-^{xiti, X2t2, xsts, X4ti) is given only by graph 
(a). We have not found a complete proof that it is independent 
of external times, but we can prove the weaker 

Lemma 1: 



^uiuui^iti, X2t2, X2.h, Xiti) (4.8) 



(see graph (a) on figure 2) is independent of the most advanced 
time. 

Proof; 

First suppose that a response-function enters at the most- 
advanced time t' . Then there is the following cancellation 




= (4.9) 



The mechanism is the same as in the proof of theorem 1 ; since 
by assumption t' is the most advanced time, the argument of 
the right-most disorder can never change sign, and can be in- 
tegrated over, even though it is odd, i.e. A'(0+). Thus the 
remaining diagrams have the structure 




(4.10) 



This diagram is independent of t', as long as t' is the largest 
(external) time. 

Note that for a loop made out of two disorders, there is only 
one diagram remaining, namely 




It is manifestly time independent. 
Lemma 2; 

r^fifiC^^i^i v^2i2, x:,h,XiU) (4.12) 

(see graph (a) on figure 2) is independent of the differences in 

time, if those are very large. 

Proof: 

By inspection, one finds that by taking the external times, 
i.e. one time at each disorder, (infinitely) far apart, the remain- 
ing integrals become unambiguous. Thus integrating over the 
response-functions does not leave any time-dependence. 
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<II1> 




FIG. 4: The two contributions to {u^) ^ at leading order. 



As mentioned above, calculation of the diagram (a) be- 
comes then possible and one finds the 

Property (missing acausal loop): 

The diagram (a) is given by (—1) times the acausal loop, 
if there one replaces each response-function by a correlation- 
function. 

Intuitively this means, that if the acausal loop would give a 
contribution, then all diagrams would cancel. This seems to 
be a general mounting theorem. 

Check: 

One can calculate diagram A on Fig. 4 explicitly using 
none of the above theorems or conjectures. The result is a 
formidable expression, which has to be integrated over mo- 
menta. By evaluating it for given values of the momenta (not 
even necessarily conserving momentum), one can compare 
with the prediction of property 5. We found both expressions 
to be equal for any randomly chosen values of the momenta. 

The properties described here suggest the following 

Property (any loop): 

All graphs can be computed to any number of loops, using 
generalizations of the above rules. 

This is not attempted here but prehminary investigation 
suggests that the same mechanism holds for two loops with 
some simple end result related to the signs of possible 
"fermion loops". 

V. FINAL RESULT FOR THE FOURTH 
CUMULANT AND UNIVERSAL RATIO 

In this Section we compute the fourth cumulant 



D 



J XV 



as well as the ratio (kurtosis): 



R 



D 



2 f 

Jxy 



(5.1) 



(5.2) 



which, according to the discussion in Sections II and III is 
expected to be universal and characterizes the deviations from 
the Gaussian approximation (for which R = Q). 

Below, we compute R at depinning both for short range (re- 
spectively long range) elasticity to lowest non-trivial order in 
e = 4 — d (respectively e = 2 — cf), i.e within a 1-loop cal- 
culation. However, since it turns out that the momentum inte- 
grals involved in the calculation depend very strongly on the 
dimension, we found it useful, and sometimes more accurate. 



to carry a 1-loop approximation directly in fixed dimension d. 
Also, since there is one exact result for a massive propagator, 
we also give the result in that case. 
We denote: 



9{D = ^ 



(5.3) 



with the obvious change for long range elasticity g{q) = , 
and (see below) the massive propagator. 

The final result in the continuum is given by the sum of the 

two diagrams in Fig. 4: 



Jxy 



d'^q d'^k d'^p 

MQfgikfgipfgip + <i)9{p + k) 

+g{qfg{kfg{p)g(p + q)g{p + k)g{p + k + q) 



(5.4) 



d'^q_ 



(5.5) 

The combinatorics can be done as follows. There is a factor 
l/(4!2^). There are 4! ways to associate each one of the four 
external u to an unsplitted vertex. Say 1, 2 are now linked to 
and 3, 4 to M^. At each vertex a u field comes out. There 
are 2* ways to assign them to each splitted vertex. Then there 
is a unique set of four splitted points, one at each vertex, en- 
tering the acausal graph (which - in effect - is the only one 
arising, with the minus sign). But there are still 3 ways to join 
these four points in a loop: Two give the first integral, one the 
second, and finally for each case the orientation can be chosen 
in two ways. 

Let us go to the discrete model with periodic BC [34]. We 
recall that 

where here and in the following the term with n = is always 
excluded. In the limit of large L/a one finds: 



(5.7) 



-2A'(0+)^i-^''( 



E 

+ 



\16 



2tt' 



1 



1 



(n2)2(TO2)2(p)(; + n)2(; + m)2(/ + m + n)^ 



2n 



y- 



(5.8) 
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One can see that the ratio R will be universal since the 1- 
loop FRG fixed point equation taken at u = yields: 



(e - 2C)A(0) = (e/)A'(0+)2 
(eJ) = LOlI ^ 



(47r)'^/2r(3- 87r2 



(5.9) 
(5.10) 



The last identity is valid for d = 4. In fact, since the 1-loop 
FRG equation is universal, it holds as well for d = 4 as for 

d < 4. For d < 4, we use 

f d^p 1 1 
U ^ 1 



LdrJ = el = e 



y- 



(5.11) 
(5.12) 



This is all we need to compute the universal ratio. 



A. Calculation to lowest order in e = 4 - d 



The ratio i? = D/(2E,,G2^) is: 



E 



1 



1 



(27r)^ J2neZ'' n» 

1 



1 



(n2)2(TO2)2(Z2)(; + nf^i + rnf{i + ^ + nf _ 

(5.13) 



One finds, using 



^e-*"' =6(3,0, e- 



(5.14) 



that in d = 4 



1 1 

^ ^ = g y i'(e(3, 0, e-*)'^ - 1) = 10.2454 (5.15) 



Noting /(p) = Y.a^7d in2yiL^„yi One has: 



E 



E 



1 



E (^/W'-1850. 



(5.16) 



(n2)2(m2)2(Z2)(/ + nf{l + mf{l + m + nf 
w 980. (5.17) 
The final result is: 



R = -l.n^e^ « -0.13e2 



(5.18) 



This results shows that R is quite small near d = A, but in- 
creases quite fast as the dimension is lowered. However the 
sums over the momenta depend very strongly on d (see below) 
and one should expect significant higher order corrections in 
e. Thus the result (5.18) is hkely to drastically overestimate 
the (absolute value of the) result in lower dimensions, which 
is why we now turn to an estimate in fixed dimension. 

B. 1-loop estimate in general dimension 

One can perform an estimate in general dimension, based on 
an arbitrary truncation on (i) 1-loop graphs (ii) neglect of the 
finite size scaUng function. 

In general dimension one has the dimensionless ratio: 



R = -{e-2Cf 

E 



1 



2 

(n2)2(m2)2(i2)2(; _^ „)2(; + 

1 



(n2)2(TO2)2(;2)(; + nf{i + to)2(2 + ^ + 

(5.19) 



Let us give a table of values: 



~ 45 



2.16465 (d = 1) 



= 6.0268120 
= 16.5 



y- 



E 



4725 



2.00815 



4.28143066080578 
6.9458079272 



(d = 2) 
(d = 3) 



(d=l) 

(d = 2) 
(d = 3) 



(5.20) 



(5.21) 



(n2)2(m2)2(Z2)2(; + nY^i + „)2 

= 0.37342751117 (d = 

= 26.567 (d = 

= 240 (d = 



1) 
2) 
3) 

(5.22) 



E 

n,m.,lGl 



(n2)2(m2)2(/2)(/ _l_ „)2(; _|_ „i)2(-/ -\-rn + ny 



= 0.069672062794960 
= 14.138 
= 143 



In d = 1 one finds: 



R = -0.0867759691287 (^1-2^ 

« -0.00964177 (C = e/3) 
« -0.00347 (C = 1.2) 



(d=l) 
(d = 2) 
(d = 3) 



(5.23) 



(5.24) 

(5.25) 
(5.26) 
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In = 2 one finds: 



R = -0.4326 I 1 - 2^ 



-0.0481 (C = e/3) 



(5.27) 



In d = 3 one finds: 



R = -0.3297 1 -2 



X 



-0.0366 (C = e/3) 



(5.28) 



where we have inserted various choices for C, including the 
1-loop result C = e/3. One finds that already in d = 3 the 
1-loop approximation is significantly lower than the extrapo- 
lation from (5.18) as discussed above. 

C. Long range elasticity in general dimension 

For LR elasticity the upper critical dimension is d^c = 2. The 
general expression for R is: 



R = -{e-2Cf 

E 



2 



{n^){m^){F)\l + n\\l + m\ 



1 



{n'^){m'^)\l\\l + n\\l + m\\l + m + n\_ 

(5.29) 



It is interesting to compute Rmd= 1. Using that: 

(5.30) 
(5.31) 

3.847 (5.32) 



V ^ = — = 3.28987 
^ n2 3 

Y — = ^ = 2.16465 



E 



1 



E 

n,m.,l^Z 



one finds: 



{n^){m'^){P)\l + n\\l + m\ 
1 

{n'^){m'^)\l\\l + n\\l + m\\l + m + n\ 



1.934 



(5.33) 



R = -0.4109 1-2 



w -0.04566 (C = e/3) 
D. Long range epsiion expansion 



(5.34) 



Similarly, one can perform an expansion in e = 2 — d. In 
d = 2 one has: 



e/= l/(27r) 



(5.35) 



45 



El tt' 

E 



: 6.0268120 



{n'^){m'^){P)\l + n\\l + m\ 



E 



n,m,leZ^ 



(n2)(m2)|/||Z + n\\l + m\\l + m + n 



(5.36) 

550 ± 20 (5.37) 

= 370 ± 10 
(5.38) 



This yields 



i? = -e^ (1 - 2Ci)' (27r)2 ^ 



E 



1 



+ 



{n^){m,^){l2)\l + n\\l + m\ 
1 



(n^)(m^ 



n||l + m||l + m + n| 



The result is 



R 



-6.17-e2 
9 



-0.686e^ 



(5.39) 



(5.40) 



again, a probable overestimation of the result if naively ex- 
trapolated tod= 1. 

E. Harmonic well, SR elasticity 

It is interesting to compare with the calculation in a massive 
scheme, i.e. an interface in a harmonic well. Setting g{q) := 
1/(1-1- g'^), we have for d = 1 



g{qf = - = 1.5708 
g{qf = ^ = 1.62596 



/// 

J q J h J p 



16 

9{<lf9{kf9{pfg{p + q)g{p + k) 



163l7r3 



= 1.62588 



(5.41) 
(5.42) 

(5.43) 



/// 

J q J k J p 

_ 2457r3 
~ 5184 
This gives the ratio 
2366 



31104 

g{<lf'9{kf9{p)9{p + Q)g{p + k)g{p + k + q) 



1.46538 



(5.44) 



R = — 



(1 - 2C/e)2 = = -0.216369455 



1215' " ' 10935 

(5.45) 

It is interesting to compare the present result to one case 
(to our knowledge the only one apart from mean field models) 
where the full distribution ofu is known in a non- trivial disor- 
dered problem [30]. This is the static random field model in 
(i = in a harmonic well (i.e. the massive case), the so-called 
toy model. The exact result there is: 



= -0.080865 . . . 



(5.46) 
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The present 1-loop approximation for the problem of de- 
pinning, continued to d = would give the larger result 
R = —1/3. It is unclear at present whether the difference be- 
tween the two results indicates that the 1-loop approximation 
is unsatisfactory so far from d = 4, or if statics and depinning 
have radically different values of R. 

VI. BEHAVIOR AT THE CRITICAL 
DIMENSION 

In this Section we reexamine isotropic depinning, and stat- 
ics, exactly in d = 4. We solve the RG equations in d = 4 
and obtain the behavior of the correlation function. Contrarily 
to periodic systems at the upper critical dimension [31], non- 
periodic objects such as interfaces submitted to either random 
bond or random field type disorder exhibit a roughness, which 
is a power of a logarithm. 

The FRG flow equation (/3-function) for the (renormalized) 
force correlator has a good limit for d = 4. If one defines 



Aiiu) = 87r2;2Ci-iA,(u/-Ci) 



(6.1) 



with I = ln(A/m) (A some UV cutoff) then the function 
A;(u) satisfies 



diA(u) = (1 - 2Ci)A(m) + CiuA'iu) 
1 

~2 



+ 1-^I32{A) . (6.2) 



(A(«)-A(0))^ 
where, for depinning, 

/32(A) = [(A(,i) - A(0))A'(w)2]" + A'(0+)2A"(,i) . 

and Ci = C/c — 1/3 is the 1-loop value, see e.g. [10]. It is 
then easy to see that the function A; {u) converges towards the 
1-loop fixed point with the foUowing asymptotic corrections: 

Ai{u) = A*{u)+J2l~^"Uu) + j(}[[AT'p2iA*) 

n 

(6.3) 

where (J3[)^^ is the inverse of the linearized 1-loop (3- 
function and the a>„ are the 1-loop eigenvalues. 

Using (6.1) with = 1/3 yields the result for the correla- 
tion function at g' = 0: 

{ugU_g)\g^rn = m-^A(O) (6.4) 
= cm"^ ln(A/m)-3 (1 + 0(1 / ln(A/m))) 

as m ^ 0, with c = 87r^A*(0), both for statics and depin- 
ning; the difference hes in the subdominant piece. Within the 
present approach using the renormaUzation scheme atq = 0, 
the 2-point correlation function at non-zero q can be computed 
from the renormalized (uniform) effective action by resum- 
ming an infinite set of diagrams. Using the standard finite size 
scaling ansatz allows to obtain the other Umit of the scaling 
function, where A » g » m. To lowest order (one loop) in 
the renormalized disorder one has [11] 

(g2 + m^fiugu.g) = (A(0) - A'{0+f{I{q) - 7(0)) + . . .) 

^^^^ (p'+m2)((p + (7)2+m2) • ^^-^^ 



This gives: 

{ugU-g} = ciq"^ + m^y^ ln( — ] |1 




3ln(^) 



+ ... 



Assuming scaling, i.e. that the function (1 — + . . .) 
{1 + x)~i one finds that: 



(6.6) 



and thus: 



{Ugu.g)^q-' [in(^) 
(u,-uo)2~(lna;)2/3 



1/3 



(6.7) 



(6.8) 



VII. STABILITY OF THE 1 -LOOP FIXED 
POINT 

Here we analyze the stabiUty of a functional fixed point. Two 
cases have to be distinguished: 

(a) There is the freedom to rescale the field u while at 
the same time rescaling the disorder correlator. This 
includes the random bond and random field interface 
models. 

(b) There is no such freedom, since the period is fixed by 
the microscopic disorder This is the case for a charge 
density wave (random periodic problem), but also for 
the random field bulk problem, in its treatment via a 
non-Unear sigma model. 

We first analyze the simpler case (b). 

A. Stability of a functional fixed point; periodic 
case 

Be the flow-equation given by 

diR{u) = p[R]{u) = eR{u) + f[R, R]{u) , (7.1) 

where / is some biUnear form of R, which contains at least 
one derivative for each R. A similar equation of course exists 
for A(m) = —R"{u) and the corresponding /3[A](u). 

Suppose R* (u) is the non-trivial fixed point of order e, i.e. 
P[R*] = 0. Two eigenfunctions and eigenvalues above R* {u) 
can be identified. 

(i) The constant mode dR{u) = 1 with eigen-value e. (As 
long as it is permissible physically.) 

(ii) The first subleading eigenfunction 5R{u) = R*{u) 
with eigenvalue — e. 

Proof: 



For case (i), we have for k <^ 1 

di{R*{u) + k) = ^[R* + k]{u) = eK , (7.2) 

since / does not couple to the constant by assumption. This 
proves (i). 
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For case (ii), we have {k <C 1): 

di {R*{u) + kR*{u)) = (3[R*{1 + k)]{u) 

= eR*{u){l + k) + (1 + Kff[R*, R*]{u) (7.3) 

Subtracting /3[i?*](u) = on the r.h.s. and expanding for 
small K yields 

di {R*{u) + kR*{u)) = eKR*{u) + 2Kf[R*,R*]{u) 

= -eKR*{u) , (7.4) 

where in the last equation we have again used the fixed-point 
condition l3[R*]{u) = 0, i.e. eR*{u) + f[R*,R*]{u) = 0, to 
eliminate f[R*, R*]{u). This proves (ii). 

Using the same line of arguments, it is easy to see that when 
starting from R{u) ^ R* (u) with some arbitrary ampUtude, 
the flow is always remaining on the critical manifold spanned 
by R*{u). 

Of course, there are in general more eigen-functions and 
eigen-values. See [10] for an exphcit example. 

B. Perturbations of the fixed point in presence of 
the freedom to rescale 

We state the following 
Theorem: 

The differential equation of the form 

-mdmAiu) = /3[A] 

/3[A] = (e - 20 A(«) + CuA'iu) + /[A, A] (7.5) 

where the symmetric functional /[A, A] transforms under 

A(u) — > K^^A(K?i) in the same way as A, has the two eigen- 
functions and eigenvalues of perturbations around the fixed 



point /3[A*] = 

zoiu) = wA'(w) - 2A(u) (7.6) 

Ao = . (7.7) 

zi{u) = CuA'{u) + {e- 2C)A(m) (7.8) 

Ai = -e . (7.9) 



(We note A instead of A* for the fixed point for simplicity 
of notations.) Note that the assumptions are satisfied by the 
1-loop flow-equation (RF-case). 

Proof: 

Be (3[A]{u) = 0. Due to the assumptions, for all k 

P[k-^A]{ku) = . (7.10) 

Deriving w.r.t. k gives with (7.5) at k = 1 

(e - 20 (uA'(u) - 2A(u)) + (u (uA'(u) - 2A{u))' 
+2f {A{u), uA'{u) - 2A{u)) = (7.11) 

This is nothing but the eigenvalue equation for the perturba- 
tion zo{u) about the fixed point /3[A] = and proves the 
existence of the solution zo{u) with eigenvalue Aq = 0. Note 
that this redundant operator with eigenvalue persists to all 
orders in perturbation theory. 



Let us turn to the next solution. Multiplying (7.11) with C, 
and adding 2e times /3[A] =0 gives 

(e - 20 (CuA'(?i) + (2e - 20A(w)) 
+ Cu (CuA'(w) + (2e - 20 A(u))' 
+ 2/ {A{u)XuA'{u) + (e - 20A{u)) = , (7.12) 

where we used the bilinearity of A{u). Rearranging yields 

(e-20 {CuA'{u) + (e-20A{u)) 
+ Cu{C,uA'{u) + {€-20A{u))' 
+ 2/ (A(m), CuA'iu) + (e - 2QA{u)) 
= -e(CuA'(u) + (e - 20A(w)) . (7.13) 

This equation is nothing but the eigenvalue-equation for the 
perturbation zi (m) = (uA'{u) + (e—2C,)A{u) about the fixed 
point P[A] = with eigen-value Ai = — e. q.e.d. 
Remark: 

Consider the case of short-range disorder, i.e. that A{u) falls 
off rapidly, and is monotonic. Usually, the leading fixed point 
solution has no knot (no u such that A{u) = 0). Then zq{u) 
has no knot and zi (u) has one knot. Eigenvalues should be 
ordered (like in quantum mechanics) due to their number of 
knots. Thus we should have found the two leading solutions 
for short-range disorder. This is confirmed by the numerical 
analysis given in subsection Vll C. 

Corollary: (Random bond case) 
The differential equation of the form 

-mdmR{u) = [3[R]{u) (7.14) 
(3[R]{u) = (e - 40i?(u) + C,uR'{u) + f[R,R]{u) 

where the synnmetric functional f[R,R] transforms under 
R{u) k~'^R{ku) in the same way as R, has the two eigen- 
functions and eigenvalues of perturbations around the fixed 
point f3[R*] = 

zoiu) = uR'{u) ~ AR{u) (7.15) 

Ao = . (7.16) 

zi{u) = CuR'{u) + {€- 40i?(u) (7.17) 

Ai = -e . (7.18) 

Note that the assumptions are satisfied by the 1-loop flow- 
equation (random bond case). 

Proof: 

This can either be proven along the same lines as for (7.5) 
or by deriving zq (u) and zi (u) twice w.r.t. u and then using 
the theorem for the random field case (7.5). 

C. Numerical analysis of the RF-fixed point 

We start from the 1-loop flow equation 

-mdmA{u) = (e - 20A(u) -|- CuA'{u) 

-l[{A{u)-AiO))^]" . (7.19) 
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1, 
0.8- 
0. 6- 
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-0.2 



FIG. 5: The solution zi(u) (in black, with 2:(0) — 1) and asymptotic 
behavior as given in (7.24) (blue/bright). 



It has the following solution [5, 10] 



yi{u) - Inyi(u) = 1 + 



(7.20) 



Perturbations around this solution satisfy the differential equa- 
tion 



-mdm [A(u) 
(l-3A)z(u) 

= 



~ \ez{u) 
uz'{u) - [{yi{u) 



(7.21) 

2/i(0))(z(u)-z(0))]" 
(7.22) 



In order to have a criterion for the numerical integration, one 
has to determine the behavior at infinity. Using that for u 

oo 



yiiu) 



(7.23) 



and assuming exponentially fast decay for z{u), one finds that 
the asymptotic behavior is 



z{u) w z(0) 



(^2 - i)e-i-"V2 
2 + 3A 



(7.24) 



Quite surprisingly, the asymptotic behavior is fixed, including 
its amplitude[38]. In any case, slower power law decay is 
ruled out on physical ground since we are considering short 
range disorder. 

The solutions for A = and A = — e are given in Eqs. (7.6) 
ff. The solution zi [u) is dominant, and gives the correction 
to scaling exponent w = — e. Note that this exponent is mi- 
nus the engineering dimension of the coupling, as is the case 
in standard field-theory[32], and also for the random periodic 
class [10]. Ref. [26] cites the value oj = — e/3. We find the 
corresponding numerical solution to decay as u^^, incompat- 
ible with (7.24) and physically unacceptable [34]. 

The question arises, of whether there are more solutions 
with fast decay. Intuitively, one would expect this: Making 
A more negative, the solution overshoots and one might think 



of fine-tuning A such that it approaches the axis for large u 
from above. However this is incompatible with the asymptotic 
behavior in Eq. (7.24), which predicts, that all solutions for 
A < —2/3 converge from below. In fact, we have not been 
able to find any further solution, and we conjecture that there 
is none. It would be interesting to prove this rigorously. This 
behavior is in contrast to the random periodic case, solved in 
[10], which has infinite many subleading contributions. 

VIII. CONCLUSION 

In this paper we have explored further properties of the field 
theory of depinning. We have defined and computed universal 
observables, such as the distribution of the interface width and 
the ratio R of the connected 4-point cumulant to the square of 
the 2-point one (kurtosis). This ratio measures the deviations 
from a Gaussian approximation which we have also used to 
obtain the universal distribution. Higher order connected cu- 
mulants can be computed in a similar way to one loop using 
polygon diagrams, and one should be able to reconstruct the 
full distribution from them. Other properties of the theory 
such as the behavior at the upper critical dimension and the 
finite size scaling behavior have been clarified. All calcula- 
tions in the present paper are of interest for comparison with 
numerical simulations, existing ones [16] or in the near future. 

In the process of computing the 4-point function we dis- 
covered massive cancellations between diagrams. We traced 
this back to the physically expected property that correlations 
exactly at depinning should be time independent. A diagram- 
matic proof of this property is still incomplete, but we have 
provided some convincing elements in that directions. As a 
result the correlation functions can be computed in a much 
simpler way. Thus there seems to be an underlying theory, 
with "quasi-static" diagrams (i.e. not containing time explic- 
itly), with some additional rules. We have understood these 
rules to lowest ( 1 -loop) order and it would be of high interest 
to understand -and prove- them to all orders. It is even pos- 
sible that there exists a simpler formulation of the theory at 
depinning in terms of, e.g. effective fermions. The fermionic 
character is suggested by the cancellation of all diagrams ex- 
cept for the "acausal loops" with a minus sign. We thus en- 
courage further examination of this fascinating question and 
full elucidation of the field theory which describes depinning. 

APPENDIX A: DIFFERENT BOUNDARY 
CONDITIONS 

a. Periodic boundary conditions 

A periodic function u with period L satisfies 

u{x + Lei) — u{x) , (A.l) 
where the ii are orthonormal. It can be written as 



Uk e ' 



(A.l) 



where summation runs over all k, such that ki = niiii/L, 
rii g Z. 
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b. Open boundary conditions 

To simplify the notation, we give all formulas for one dimen- 
sion; generalizations are straightforward. 

Suppose the function f{x) is defined on [0, L]. Then a func- 
tion g{x) can be defined by the following prescription 

/ \ _ / /(a;) for 0<x<L 

9{x) - I ^^2L -x) for L<x<2L 

g{x) can be prolonged to a periodic function with period 2L, 
i.e. g{x+2L) = g{x), and satisfies by construction in addition 



g{2L -x)= g{x) . 



(A.4) 



In the basis needed to construct functions with period 2L, we 
have to restrain ourselves to 

since the sin do not satisfy (A.4). The such constructed set 
of functions f{x), x G [0, L] has Neumann-boundary condi- 
tions at X = and x = L. From (A.5), we infer that the 
number of modes is reduced by a factor of 2 (compared to 
the case of closed boundary conditions), but the construction 



does not change any observable constructed from w"^ or any 
energy, all based on the symmetry relation (A.4). Importantly, 
the modes have all mean 0, which is not the case for other 
basis, e.g. when using anti-periodic functions. Also note, that 
this ansatz reproduces the formula in [33]. As an interesting 
consequence, we observe that the following systems lead to 
the same distribution 

• An elastic object with TV degrees of freedom, and closed 
boundary conditions. 

• An elastic object with 2N degrees of freedom, and open 
boundary conditions. 

The simplest example is a 1 -dimensional random walk with 
closed boundary conditions, and a 2-dimensional random 
walk with open boundary conditions, which can be checked 
numerically [39]. 
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